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Abstract. A theory of spin-transport in hybrid normal metal - ferromagnetic electronic circuits is de- 
veloped, taking into account non-collinear spin-accumulation. Spin-transport through resistive elements 
is described by 4 conductance parameters. Microscopic expression for these conductances are derived in 
terms of scattering matrices and explicitly calculated for simple models. The circuit theory is applied to 
2-terminaI and 3-terminal devices attached to ferromagnetic reservoirs. 
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1 Introduction attention. Moodera et al. measured the dependence of the 

current through magnetic tunnel junctions on the rela- 
Spin-injection from a ferromagnetic metal into a non-magnetifeive angle between the magnetization directions of the 



metal was first realized by Tedrow and Meservey in the 
seventies The discovery of the giant magneto-resistance 
(GMR) in metallic magnetic multilayers has led to an ex- 
plosion of interest into spin-transport in hybrid normal 
metal-ferromagnetic metal systems in the nineties (for re- 
views see Ref. Q). Although discovered only ten years 
ago, the GMR is commercially utilized in high-end mag- 
netic recording media. Spin-transport between ferromag- 
nets through tunnel junctions has also attracted renewed 
interest (for a review see Ref. ^). Recently a magnetic 
double barrier tunnel device was fabricated which enables 
the study of the interplay between spin-polarized tunnel- 
ing and Coulomb charging effects^,^,||. 

Magneto-electronic principles have until now led to 
magnetic-field sensor devices. Future applications might 
include non-volatile memory cells or even transistors, the 
latter being a three-terminal device. Johnson realized that 
novel long-range effects in transport between ferromagnets 
and normal metals can occur in multi-terminal systems . 
In Johnson's spin-transistor different ferromagnetic (Ohmic) 
contacts provide information about the amount of spin- 
accumulation in a normal metal film over distances much 
larger than the mean-free path. A transistor-like effect was 
observed on switching the configuration of the magnetiza- 
tions of the ferromagnetic contacts from parallel to anti- 
parallel. 

Spin-transport in systems comprising ferromagnets with 
non-collinear magnetization directions has attracted less 
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electrodes [g. The results were in agreement with the the- 
oretical predictions based on the spin-torque model intro- 
duced by Slonczewski|^. Non-collinear spin transport has 
been addressed by a number of theoretical papers^]. 

The above examples do not exhaust the novel phe- 
nomena that can be anticipated for multi-terminal hybrid 
normal metal-ferromagnet circuits and devices. The ques- 
tion of the most appropriate theoretical approach to the 
field arises. If the possibilities of macroscopic quantum 
coherence and its application to quantum computing are 
contemplated, a fully quantum mechanical treatment of 
the many-body system is required, of course. However, 
for contact with the experiments mentioned above and 
probably also with most to be realized in the near future 
a full quantum mechanical treatment is unnecessary and 
unrealistic. When at least a part of the device is diffu- 
sive simplified approaches are called for. The situation is 
similar to that in the field of inhomogeneous supercon- 
ductors which has recently been reviewed by Belzig et 
al. |l^. The theoretical framework of choice for trans- 
port in superconducting and/or magnetic "dirty" systems 
is the non-equilibrium Keldysh Green functions formalism 
in the quasiclassical approximation. However, it is techni- 
cally difficult and physically not very transparent for all 
but the devoted specialist. This led one of us to simplify 
equations for complicated hybrid superconductor-normal 
metal systems to a handful of easily accessible rules, the 
circuit theory of Andreev reflection [^. We recently in- 
troduced a circuit (or finite-element) theory of spin and 
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charge transport in hybrid ferromagnet- normal metal sys- 
tems which provides parametric dependencies of the elec- 
tron transport properties as well as microscopic expres- 
sions for the parameters and illustrated its appeal |p^. 
Fig. shows an example of a typical many-terminal con- 
figuration of present interest. The main physical parame- 
ters of such a circuit are the magnetization directions, the 
chemical potentials of the leads, and the contact conduc- 
tances between the normal and ferromagnetic metals. In 
this paper we introduce the spin-polarized kinetic equa- 
tions based on the Keldysh Green function technique in 
the quasi-classical approximation. When solved with judi- 
ciously chosen boundary conditions the basic equations of 
the circuit theory emerge. 

Before turning to the technical details let us first dis- 
cuss the conditions under which long-range spin effects are 
observable in normal metals. Spins injected into a normal 
metal node relax due to unavoidable spin-flip processes. 
Naturally the dwell time on the node must be shorter 
than the spin-flip relaxation time in order to observe non- 
locality in the electron transport. For a simple ferromagnet 
(F) normal metal (N) double heterostructure (F — N — F) 
with anti parallel magnetizations the condition can be 
quantified following Ref. The spin-current into the 
normal metal node is roughly proportional to the particle 
current, e{ds/dt)tr ^ I — V/R, where s is the number 
of excess spins on the normal metal node, V is the volt- 
age difference between the two reservoirs coupled to the 
normal metal node, and R is the F — N contact resistance. 
When the node is smaller than the spin-diffusion length, 
the spin-relaxation rate is e{ds/dt)rei = ^•^/''sf) where 
Tgj denotes the spin-relaxation time on the node. (Other- 
wise this simple approach breaks down since the spatial 
dependence of the spin-distribution in the normal metal 
should be taken into account p3| ) . The number of spins on 
the normal metal node is equivalent to a non-equilibrium 
chemical potential difference Z\/i = s5 in terms of the en- 
ergy level spacing 5 (the inverse density of states) (more 
generally the relation between and s is determined 
by the spin-susceptibility ) . The spin-accumulation 
on the normal metal node significantly affects the trans- 
port properties when the non-equilibrium chemical po- 
tential difference is of the same order of magnitude or 
larger than the applied source-drain voltage, A^j. > eV or 
6Tgf/h > R/Rk, where Rk = e^/h is the quantum resis- 
tance. We see that spin-accumulation is only relevant for 
sufficiently small normal metal nodes and/or sufficiently 
long spin-accumulation times and/or good contact con- 
ductances. The implications of this relation have been dis- 
cussed in Ref. js) for different materials. 

In the present article we explain in more detail the 
foundations of the circuit theory of spin-transport |^ 
and present more applications; a general result for two- 
terminal devices and results for the spin-resistance in three- 
terminal devices similar to the set-up of Johnson Q]. The 
manuscript is organized in the following way. In Section ^ 
we describe the basic entities in a circuit theory, the nodes, 
the contacts and the reservoirs. The contact conductances 
are computed in section |^ for a diffusive, a ballistic and 



a tunnel contact. The circuit theory is employed in sec- 
tion ^ to find the current trough two-terminal systems 
and the 'spin-resistance' of a three terminal device. Our 
conclusions can be found in section Appendix ^ gives 
the detailed derivation of the transformation from a non- 
collinear to a coUinear two-terminal system. 



2 Circuit theory 

A typical (magneto)electronic circuit as schematically shown 
in Fig. |l| can be divided into contacts (resistive elements), 
nodes (low impedance interconnectors) and reservoirs (volt- 
age sources). The present theory is applicable when the 
contacts limit the electric current and the nodes are char- 
acterized by a distribution function which is constant in 
position and isotropic in momentum space. The latter con- 
dition justifies a diffusion approximation and requires that 
the nodes are either irregular in shape or contain a suf- 
ficient number of randomly distributed scatterers. State- 
of-the art magneto-electronic devices are rather dirty, so 
this does not appear to be a major restriction for the ap- 
plicability of the theory. Because the spin-accumulation is 
not necessarily parallel to the spin-quantization axis, the 
electron distribution at each node is described by a 2 x 2 
matrix in spin-space. The current through each contact 
can be calculated as a function of the distribution matri- 
ces on the adjacent nodes. The spin-current conservation 
law then allows computation of the circuit properties as a 
function of the applied voltages. The recipe for calculating 
the current-voltage characteristics can be summarized as: 

— Divide the circuit into nodes, contacts, and reservoirs. 

— Specify the 2x2 distribution matrix in spin-space for 
each node and reservoir. 

— Compute the current through a contact and the dis- 
tribution matrices in the adjacent nodes, which are 
related by the spin-charge conductances specified be- 
low. 

— Make use of the spin-current conservation law at each 
node, stating that the difference between in and out- 
going spin-currents equals the spin-relaxation rate. 

— Solve the resulting system of linear equations to ob- 
tain all currents as a function of the reservoir chemical 
potentials. 



2.1 Node 

We denote the 2x2 distribution matrix at a given energy e 
in the node by /(e), where hat denotes a 2 x 2 matrix in 
spin-space,. The external reservoirs are assumed to be in 
local equilibrium so that the distribution matrix is diag- 
onal in spin-space and attains its local equilibrium value 
/ = lf{e,^a), 1 is the unit matrix, f{e,iJa) is the Fermi- 
Dirac distribution function and /Zq is the local chemical 
potential in reservoir a. The direction of the magnetiza- 
tion of the ferromagnetic nodes is denoted by the unit 
vector nia. 
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The 2x2 non-equilibrium distribution matrices in the 
nodes in the stationary state are uniquely determined by 
current conservation 



processes in the contacts, the reflection matrix for an in- 
coming electron from the normal metal transforms as 




where lap denotes the 2x2 current in spin-space from 
node (or reservoir) a to node (or reservoir) (3 and the 
term on the right hand side describes spin-relaxation in 
the normal node. The right hand side of Eq. (|^) can be 
set to zero when the spin-current in the node is conserved, 
i.e. when an electron spends much less time on the node 
than the spin-flip relaxation time Tgj-. If the size of the 
node in the transport direction is smaller than the spin- 
flip diffusion length l^^ — yHDr^, where D is the diffusion 
coefficient then the spin-relaxation in the node can be in- 
troduced as {df^/dt)j.^i = (iTr(/^)/2 - /^)/rgf. If the 
size of the node in the transport direction is larger than 
Zgf the simplest circuit theory fails and we have to use a 
more complicated description with a spatially dependent 
spin-distribution function|fl3[. 



2.2 Current through a contact 

A schematic picture of a contact between a normal metal 
and a ferromagnetic node is shown in Fig. |^. The current is 
evaluated on the normal side of the contact (dotted line). 
The current through the contact is 



(2) 



where r"™ (s —1) and r""' (s =|) are the spin-dependent 
reflection coefficients in the basis where the spin-quantization 
axis is parallel to the magnetization in the ferromagnet, 
the spin-projection matrices are 



(i + 6--m)/2, 



(i 



m)/2 



(3) 
(4) 



and is a vector of Pauli matrices. Similarly for the trans- 
mission matrix 



where t""^ and t"™ are the spin-dependent transmission 
coefficients in the basis where the spin-quantization axis 
is parallel to the magnetization in the ferromagnet. Using 
the unitarity of the scattering matrix, we find that the 
general form of the relation (p|) reads 



el - G^u^ f -f'']u^+ G^u^ if 



-G^^u^fu^ - {G^^y-u^fu^ 



(5) 



where we have introduced the spin-dependent conductances 
and 



where r"™ is the reflection coefficient for electrons from 
transverse mode m with spin s' incoming from the nor- 
mal metal side reflected to transverse mode n with spin 
s on the normal metal side, and i^"™ is the transmission 
coefficient for electrons from transverse mode m with spin 
s' incoming from the ferromagnet transmitted to trans- 
verse mode n with spin s on the normal metal side. (Note 
that the hermitian conjugate in (|^) operates in the spin- 
space and the space spanned by the transverse modes, e.g. 
(f™")^^, = (f^,'^)*) The relation (|) can be found intu- 
itively in the spirit of the Landauer-Biittiker formalism, 
but can also be derived more rigorously by using the Keldysh 
formalism for non-equilibrium transport. We use below the 
latter approach and clarify the order of the spin- indices. 
Implicitly included in (||) are also effects related to the 
precession of spins non-coUinear to the magnetization di- 
rection in ferromagnets which will be made more explicit 
below. 

The relation (^) between the current and the distri- 
butions has a simple form after transforming the spin- 
quantization axis. The detailed calculation of this trans- 
formation is shown in Appendix Disregarding spin-flip 



gt = 



G^ 



M 



M 



nm 1 2 

T I ■ 



El 



^nm 1 2 



\ ^ I j.nm |2 



nm 



(6) 



(7) 



and the mixing conductance 



G 



Ti 



h 



M-^rp(rp)* 



(8) 



The precession of spins leads to an effective relaxation of 
spins non-coUinear to the local magnetization in ferromag- 
nets and consequently the distribution function is limited 
to the form = + (t • mff. Such a restriction does 
not appear in the normal metal node and can be any 
hermitian 2x2 matrix. 

We thus see how the relation between the current through 
a contact and the distributions in the ferromagnetic node 
and the normal metal node are determined by 4 parame- 
ters, the two real spin-dependent conductances (G^, G^) 
and the real and imaginary parts of the mixing conduc- 
tance G^y These contact-specific parameters can be ob- 
tained by microscopic theory or from experiments. The 
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spin-conductances and have been used in descrip- 
tions of spin-transport for a long time|2|. The mixing con- 
ductance is a new concept which is relevant for transport 
between non-collinear ferromagnets. Note that although 
the mixing conductance is a complex number the 2x2 
current in spin-space is hermitian and consequently the 
current and the spin-current in any direction given by 
Eq. (||) are real numbers. From the definitions of the spin- 
dependent conductances (^), (|^) and the 'mixing' conduc- 
tance (Isl) we find 



2ReG 



Ti - n^i 



G' +G^ 



■El 



2.3 Derivation of the current 

In this section the relation between the current through a 
contact and the adjacent distribution functions (|^) is de- 
rived. This derivation is not crucial for the understanding 
of Sections ^-|^ and can be skipped by readers mainly inter- 
ested in the physical implications of (^) and (^) . We follow 
the lines of the derivation of the contact current between a 
normal metal and a superconductor in Ref. ||l2|. We con- 
sider a F-N system comprising of a ferromagnet with ar- 
bitrary magnetization direction and a normal metal node 
separated by a contact, as shown in Fig. The Stoner 
Hamiltonian is 



and consequently the conductances should satisfy 

2ReG^^ >G^ +GK (9) 

A physical interpretation of this result is given below. 

Before outlining the derivation of (||) leading to (|^) , let 
us discuss the physics in simple terms. Some insight can 
be gained by re-writing the current and the distribution 
function in terms of a scalar particle and a vector spin- 
contribution, / 

and F = i/(f ' 
be written as 



eN 



= (l/o + a- ■ I,)/2, = If^ + <T • s/, 
(T • mff . The particle current can then 



(GT + G^)(/o''-/(f) 

+ (GT-Gi)(/f -m-s/f). 



(10) 



The familiar expressions for coUinear transport are recov- 
ered when m • s = ±1. The spin-current is 

ls=m[{G^ -G^){f^ -f^) 

+ {G^ + G^)/f + (2ReGT^ - G'f - G^)s ■ m/f ] 
-s2RcGT^/f + (s X m)2ImGT^/f . (11) 

The first three terms point in the direction of the magne- 
tization of the ferromagnet m, the fourth term is in the 
direction of the non-equilibrium spin-distribution s, and 
the last term is perpendicular to both s and m. The last 
contribution solely depends on the imaginary part of the 
mixing conductance. We can interpret this term by con- 
sidering how the direction of the spin on the normal metal 
node s would change in time keeping, all other parame- 
ters constant. The cross product creates a precession of s 
around the magnetization direction m of the ferromagnet 
similar to a classical torque while keeping the magnitude of 
the spin-accumulation constant. In contrast, the first four 
terms represent diffusion-like processes which decrease the 
magnitude of the spin-accumulation. We now understand 
condition (^) since (pH ) implies that the non-equilibrium 
spin-distribution /j^^propagates easier into a configura- 
tion parallel to s than parallel to m, since these processes 
are governed by positive diffusion-like constants 2ReG^''- 
and 2ReGT^ - G^ - G^, respectively. 



H 



1 
2m 



V2 + yP(r) 



i + l>'(r). 



(12) 



where yP(r) is the spin- independent potential and ^^(r) 
is the spin-dependent potential. The latter vanishes in the 
normal metal node and has an arbitrary, but spatially 
independent direction in spin-space in the ferromagnetic 
reservoir, y^(r) = (<t • m)l/*(r) . Non-equilibrium trans- 
port properties are most conveniently discussed in the 
framework of the Keldysh formalism. The Keldysh Green 
function is given by the (4 x 4) matrix 



G = 



/G« G^ 
I 6 G^ 



where G^, G^ and G^ are the retarded, Keldysh and 
advanced Green functions respectively, which are (2 x 2) 
matrices in spin-space and is the (2 x 2) zero matrix. 
The retarded Green function in spin-space is 



G^(l,l') 



Gf^(l,l') GUI,!') 
GfT(l,l') Gh(l,l') 



and there are analogous expressions for G^ and G^. Here 
1 denotes the spatial and the time coordinates, 1 = Titi. 
The symbol "check" Q denotes (4 x 4) matrices in Keldysh 
space and the symbol "hat" denotes (2 x 2) matrices in 
spin-space. The spin-components of the Green functions 
are 

G« (1, 1') = -^e{t, -tr){ [Mil ^l(l')] )+ , (13) 
G^,(l, 1') = ^e{t,, t,) { [Ml), ^1(1')] ) , , (14) 



Gf,(l,l') = -*([V'.(1),V'1(1')]) , 



(15) 



where V'l(l) is the electron field operator for an electron 
with spin s in the ^-direction. The Keldysh Green func- 
tions is determined by the equation 



H - ih^ ) G(rt, r't') = iS{rt ~ r't' 



and the boundary conditions to be discussed below (1 is a 
4x4 unit matrix). In the stationary situation G(l, 1') = 
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/ d{E/2n) exp{iE {ti - tr))GE{r, r') and the Green func- 
tion on a given energy shell is determined from 

(h-E) Gsivy) = \S{y-y'). 

We will in the following omit the index E in denoting 
the Green function at a given energy. The Keldysh Green 
function can be decomposed into quasi-one-dimensional 
modes as 



Gss'{vy) 



^nsms' 2; ) X 



71771. a/3 



(16) 



where x" (p; ^) is the transverse wave function and fc" de- 
notes the longitudinal wave- vector for an electron in trans- 
verse mode m with spin s. The indices a and (3 denote 
right-going (-|-) and left-going (— ) modes. The symbol 
'tilde' Q denotes matrices in Keldysh space, spin-space, 
ant the space spanned by the transverse modes and the 
directions of propagation. 

The current operator can be found from the continuity 
relation for the electron density. The spin-density matrix 
is 

= (Vl(l)^.(l)) . 
The time-evolution of the spin-density matrix reads 



prec. 



where we have inserted the Hamiltonian (|l2|) and found 
the spin-current 



and the spin-precession 



prec. 



(17) 



In three dimensions the spin-precession is an average over 
many states with different Larmor frequencies which av- 
erage out very quickly in a ferromagnet, leading to an 
efficient relaxation of the non-diagonal terms in the spin- 
density matrix that represent a spin-accumulation non- 
collinear to the magnetization in the ferromagnet. This 
spin-relaxation mechanism does not exist in normal metals 
where in the absence of spin-flip scattering the spin-wave 
functions remain coherent. 

In the stationary situation the spin-current is 



(- 


d ' 


\ hi 






dry , 


1 2rn J 





eMiEih ~ tr)) \r,,=r, ■ (18) 

We now define the extended 4x4 current matrix in Keldysh 
and spin-space as 



The transverse wave function x"(p;a;) is spatially inde- 
pendent in the leads and we find 

Ls' (x) = *e ^ (a< - /?«:?) G^ts' ^) X 

na/3 

dpx:{p;x)x:r{p;x), (19) 

where v" — fik"^ / m is the longitudinal velocity for an elec- 
tron in transverse mode n with spin s. In a normal metal, 
the transverse states and the longitudinal momentum are 
spin-independent and the Keldysh current simplifies to 



/,.Ka:) = 2ie5]az;"G^^„,,(a;,a;). 



(20) 



which we will use to calculate the spin-current on the nor- 
mal side of the contact. We use the representation 



aSa,psign{x — xi) 



(21) 



where the latter term does not contribute to the current 
on the normal side, and we have on the normal metal side 



Iss' {x) = 2e V ag""„,, (x, x) 



no. 



(22) 



We now introduce the transfer matrix M between waves 
propagating to the right (left) on the right hand side of 
the contact •f'^ {^r) and waves propagating to the right 
(left) on the left hand side of the contact i^Z) 



The elements of the transfer matrix are related to the 
reflection and transmission coefficients by 







)- 




I m ^ 


m 







Similarly we can introduce the scattering matrix 



S = 



r t' 
t r' 



so that 

S = 



— (to ) (m ) ^ 

m++ — TO^ (to )^^m ^ (m )^"'" 



where r^'^ is the reflection matrix for incoming states from 
the left in mode m and spin a to mode n with spin s, 
t^'^ is the transmission matrix for incoming states from 
the left transmitted to outgoing states to the right, r' is 
the reflection matrix for incoming states from the right 
reflected to the right, and t' is the transmission matrix 
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for incoming states from the right transmitted to the left. 
Unitarity of the 5- matrix requires S'' S = 1 and SS^ ~ 
and imphes that the transfer matrix satisfies 



(23) 
(24) 



where (17^)"^^^, = aSa^pSns.ms' is a Pauh matrix with 
respect to the direction of propagation. In order to con- 
nect the Green function to the left and to the right of 
the contact, we use the transfer matrix of the contact 

dnsms'i^ — ^2,x') = "^is" ^a" -^nsls" 9fs"1ns' — ^li^') 

and similarly for the x'-coordinate. Hence 

~g2 = Af5iMt , (25) 

where 92(1) = g{x = X2{i) , x' = 2:2(1)). Up to this point the 
Keldysh Green functions have been obtained exactly for 
the Hamiltonian (p^). Now we should include the proper 
boundary conditions to uniquely define the Green func- 
tions. To this end we introduce the assumptions of isotropiza- 
tions at the nodes. The incoming modes are assumed to 
take their quasi-classical values described by the quasi- 
classical Green functions G, whereas the outgoing modes 
are determined by the properties of the contact: |p2[ 



where gj^'^ and 2 determined by the Green function 
on the right and the scattering properties of the contact. 
The advanced Green function is related to the retarded 
Green function by 

9A = -g\i ■ 

The Keldysh component on the left side is determined 
by iM), (|^, (||), (|29|) and (M) and is dictated by the 
bounaary conditions given by the isotropization process 
to be 



(^,+.gi) (^,-Gi) =0 
{S, ~ G2) (5. + 52) = 
{S, - 32) + G2) = , 



(26) 
(27) 
(28) 
(29) 



where Gi (G2) is the isotropic Green function in reservoir 
1(2): 



(Gi 



^S"''{Gl)s 



In a normal metal, the homogeneous retarded quasi-classical 
Green function is 



Gr = 



Gr Gj^ 



1 
1 



and the advanced quasi-classical Green function is Ga = 
— Gij. (Note that here 1 means a unit matrix in the ba- 
sis of the transverse modes and spin, 1 — > 6ns.ms'i)- The 
Keldysh component of the Green function is 



G 



K.l{2) 



hi 



(2) 



1 
1 



where the 2x2 distribution matrix h is related to the (non- 
equilibrium) distribution fmictions /(e)i(2) in the nodes 

by 

hn2) = 2(2/(e)i(2) - 1) . 

The isotropization conditions (^ , (|2^) , (^) and ( ^ 
relate the retarded and advanced Green function on the 
left and right hand side of the contact: 

/ 1 0^ 



9R.1 



9R..2 



1 9b.,2 
1 



9k,i 



9 k, 2 



rhi ~gKs J 
h2r'^ h2l 



(30) 



(31) 



We now use the relation between the Green function on 
the left hand side and the right hand side of the contact 
Sf) to obtain the retarded Green functions 



= 2r 



9r,i 
9r72 



2r'. 



and the Keldysh Green function 

9k7i = t'h2t'^ + hir^ 
5++ = thit^ + r'h2r'^ . 



(32) 
(33) 



Inserting the ex pres sion the Keldysh component (j3^) into 
( p2[ ) and using (2_^ )we finally find the expression for the 
current through the contact (g). 



3 Contact conductances 

The four conductance parameters G^, G-l-, ReG^-'- and ImG^^ 
depend on the microscopic details as illustrated below for 
3 elementary model contacts: a diffusive, a ballistic, and 
a tunnel contact. 



3.1 Diffusive contact 

We consider first a diffusive contact between a normal 
metal node and a ferromagnetic node and show the rela- 
tions between the conductances (see Fig. ^). The cross- 
section of the contact is A, the length of the normal metal 
part of the contact is L^, the length of the ferromagnetic 
part of the contact is , the conductivity on the normal 
side and the spin-dependent conductivities on the fer- 
romagnetic side a^^ . The conductance of the normal part 
is G^ = Aa^ /L^ and the spin-dependent conductances 
of the ferromagnetic part are G£* = Aa^^ /L^ . The spin- 
dependent conductances of the whole contact are obtained 
simply as the diffusive ferromagnetic and normal metal re- 
gions in series: 



G 



G 



T _ 




D — 




I 




D ~ 





(34) 
(35) 
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These spin-dependent conductances {G]j and G]j) fully 
describe collinear transport (in the absence of spin-flip 
scattering). For non-collinear magnetizations the mixing 
conductance Is also needed. It can be derived from the 
scattering matrix, e.g. with the method developed in Ref. 
Here we use a much simpler approach based on the diffu- 
sion equation, describing the scattering properties of the 
contact by a spatially dependent distribution matrix. The 
current density on the normal side of the contact {x < 0) 
is i{x < 0) = <J^ dxf and consequently the total current is 

/(x<0)-G^(L^a,)/, 



Note that the component of the spin-current that is non- 
coUincar to the magnetization direction in the ferromag- 
net is not conserved on going through the interface due to 
strong relaxation induced by ([l7|). The first order differen- 
Tjtial equation ( pO| ) and the boundary conditions (p^), (pTj), 
( ^8[ ) and (|3^) uniquely determine the distribution func- 
tions and hence the conductance in the diffusive contact. 
The current on the normal side of the contact becomes 



where / is the spatially dependent distribution matrix on 
the normal side in the contact. In the normal metal node 
the boundary condition is 



(36) 



In a ferromagnet spin-up and spin-down states are inco- 
herent, and hence spins non-collinear to the magnetization 
direction relax according to (|l^) and only spins collinear 
with the magnetization will propagate sufficiently far away 
from the NF-interface. We assume that the ferromagnet is 
sufficiently strong and that the contact is longer than the 
ferromagnetic decoherence length ^ — ^D//iex, where D 
is the diffusion constant and /lex is the exchange splitting. 
The decoherence length is typically very short in ferromag- 
nets, ^=2nm in Ni wires[Q. The distribution function on 
the ferromagnetic side can then be represented by a 2- 
component distribution function 

where and it^ are the spin-projection matrices (0) and 

We allow for a spin-accumulation collinear to the mag- 
netization direction in the ferromagnet. The boundary 
condition determined by the distribution function in the 
ferromagnetic node is thus 



(37) 
(38) 



/T(. = L^) = /^T, 
f-^{x = L^) = r^. 
The total current in the ferromagnet is 

i{x > 0) = G'^^u^dxf^ + G'^^u^dxf^ . 



We assume that the resistance of the diffusive region of 
the contacts is much larger than the contact resistance be- 
tween the normal and the ferromagnetic metal. The dis- 
tribution function is in this limit continuous across the 
normal metal - ferromagnetic interface. 



/(0+) =/(o-). 



(39) 



Current conservation on the left {x < 0) and on the right 
{x > 0) of the normal mctal-ferromagnet interface dictates 



G 



N 



u\f~f'')u'+u\f~f'')u' . (41) 



The current in a diffusive contact thus takes the generic 
form (|) with = G[„ = G\j and G^^ = G^. The 
mixing conductance is thus real and only depends on the 
normal conductance. The latter results can be understood 
as a consequence of the effective spin-relaxation of spins 
non-collinear to the local magnetization direction. Those 
spins cannot propagate in the ferromagnet, and conse- 
quently the effective conductance can only depend on the 
conductance in the normal metal as ( |4l] ) explicitly demon- 
strates. 



3.2 Ballistic contact 

A simplified expression for the conductances can be found 
for a ballistic contact. Firstly, the reflection and transmis- 
sion coefficients appearing in (||) , (^ and (||) are diagonal 
in the space of the transverse channels since the trans- 
verse momentum is conserved. In a simplified model p8[| 
the transmission channels are either closed i = or open 
t — 1. The conductances (^), (0) and (||) can then be found 
by simply counting the number of propagating modes. We 
obtain the spin-dependent conductances 



G 



(42) 
(43) 



(40) 



where is the number of spin-up propagating channels 
and N-^ is the number of spin-down propagating channels. 
The mixing conductance is determined by 

G^j^ =max(Gl3,Gi,) 

and is real. In a quantum mechanical calculation the chan- 
nels just above the potential step are only partially trans- 
mitting and the channels below a potential step can have a 
finite transmission probability due to tunneling. Further- 
more, the band structure of ferromagnetic metals is usu- 
ally complicated and interband scattering exists even at 
ideal interfaces. We may therefore expect that in general 
the phase of the scattered wave will be relevant giving a 
non- vanishing imaginary part of the mixing conductance. 
First-principles calculations of the complete conductance 
matrix are therefore highly desirable. 
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3.3 Tunnel contact 

For a tunneling contact the transmission coefficients are 
exponentially small and the reflection coefficients have a 
magnitude close to one. The spin-dependent conductances 
are 



nm 



(44) 



For simple models of tunnel barriers r"™ = (5"™ exp i(f>" — 
frg™, where the phase-shift is spin- independent. We 
expand (H) in terms of the small correction (5r"™ and find 
that 

t I . 

where Gip and Gj^ are the spin-dependent tunneling con- 
ductances di^). Since the transmission coeeficients in a 
tunnel contact are all exponentially small, the imaginary 
part of G'^^^ is of the same order of magnitude as G'^^ and 
G-^^ but it is not universal and depends on the details of 
the contact. 



4 Illustrations of the theory 

We will in this chapter illustrate the appeal of the circuit 
theory of spin-transport by computing the transport prop- 
erties of two-terminal devices and a three terminal device. 
It is assumed that the normal metal node in these devices 
is smaller than the spin-diffusion length so that the spatial 
distribution function is homogeneous within the node. 



4.1 Two terminals 

First we consider a normal metal node attached to two fer- 
romagnetic reservoirs with identical contacts , e.g. G\ = 

gI = G\ Gi = Gi = G^ and G^ = gI^ = G^K The 
relative angle between the magnetization in the two fer- 
romagnetic reservoirs is 9. With the aid of (|l|) and (^ we 
find the current 



P2 tan2 6'/2 



1(0) ^-v(l 5— 

^ ^ 2 V 1 + Sgf tan^ 6'/2 -h a 



(45) 



where 



hp + 5sf (3f?K + 25sf) 
('7fl + 2<7sf)(l+5sf) 



Here, we have introduced the total conductance of one 
contact G = G^+G^, the polarization P = {G^-G^)/{G^ + 
G^), the relative mixing conductance rj — 2G^^/(G^ + 
G^) and the ratio of the 'spin- flip conductance' [|j Gg£ = 
e^/(2(5Tgf) {6 is the energy level spacing) to the conduc- 
tance of the whole device in the parallel configuration 
g(g£ = Ggf/(G/2)). The current is an even function of 
9. Note that a > f. When the magnetizations are par- 
allel (9 = 0), there is no spin- accumulation on the nor- 
mal metal node and the current is given by Ohm's law 



/p = I{9 = 0) = GV/2. The anti-parallel magnetiza- 
tion configuration {9 — tt) generates the largest spin- 
accumulation, reducing the particle current to /^p — 
I{9 = tt) = G[l - PV(1 + .9sf)]V^/2. In this case the 

magneto-resistance ratio (/p — -^Ap)/-^P ^^/(^ +5sf)' 
irrespective of the relative mixing conductance rj. Natu- 
rally, the spin-accumulation and consequently the magne- 
toresistance decreases with spin-flip relaxation time. Any 
spin flip reduces the effective polarization. The result for 
long spin-flip relaxation times (ogf ^ 1) was previously 
obtained for two tunnel junctions and we have thus gen- 
eralized it to arbitrary contacts. For general 9 the current 
depends on a and thus on the mixing conductance. In 
Fig. ^ we plot the current vs. the relative magnetization 
angle 9 for a given effective polarization P^/(l -I- g^^) and 
a number of values of a (a = 1, a = 10 and a = 100). 
In general, the current increases with increasing 77. As one 
can see from (|^) a large relative mixing conductance means 
that spins orthogonal to the magnetization in the reser- 
voirs easily can escape from the normal metal node. This 
suppresses the spin accumulation. Therefore when r] ^ 1 
and 6 is not close to tt the current approaches Ohm's law 
/ ~ GV/2. Except for the anti-parallel magnetizations, 
the angle dependence vanishes in that limit which explains 
the sharp dip at 9 close to tt. 

Fig. ^ illustrates a universal property of the two-terminal 
device with non-coUinear magnetizations that is indepen- 
dent of the contact conductances and the spin-flip relax- 
ation time. By scaling the total conductance modulation 
to the difference between the conductance in the parallel 
and anti-parallel conflguration (G/2)P^/(1 +5gf) the cur- 
rent change for any two terminal device with a diffusive 
normal metal node should be above the universal curve 
determined by the minimum value of 77, \ri\ = 1. Thus, ac- 
cording to our theory the current vs. magnetization angle 
relation for a spin valve must lie above the universal curve 
obtained for I77I = 1. 



The result ( 45 ) has been derived for normal metals is- 
lands that can be described as Fermi liquids It was re- 
cently generalized to transport through a Luttinger liquid 
by Balents and Egger|2l|. In a Luttinger liquid the current 
is non-linear in the applied source-drain bias voltage, and 
the spin-charge separation reduces the spin-accumulation[pT[ . 
Coulomb charging effects can also reduce the spin-accumulation 
in the linear response regime Q. 



(46) 4.2 Three terminals 



Let us now consider a set-up similar to the Johnson spin- 
transistor as shown in Fig. |g which was also discussed by 
Geux et al. within the context of a multi-terminal 
Landauer-Biittiker formalism for coUinear magnetization 
configurations only. 

A small normal metal node is attached to two ferro- 
magnetic reservoirs and one normal metal reservoir by 
three contacts. A voltage bias applied to the ferromagnetic 
reservoir Fl and the normal metal reservoir N causes a 
current between the same reservoirs passing through the 
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normal metal node. The spin-accumulation on the nor- 
mal metal node injected by Fl affects the chemical po- 
tential of ferromagnet F2, which is adjusted such that the 
charge current into F2 vanishes. We characterize the con- 
tact between the first (second) ferromagnet and the nor- 
mal metal node by the total conductance Gi — GI + G[ 
(G2 = + G^), the polarization Pi = (G^ - Gi)/Gi 
(P2 = {Gl-Gi)/G2) and the relative mixing conductance 
772 = 2G|V(G'J + G|) (771 = 2GP/(GI + G|)). The contact 
between the normal metal reservoir and the normal metal 
node is characterized by a single conductance parameter, 
Gsjv- is the relative angle between the magnetization of 
ferromagnet Fl and ferromagnet F2. We assume that the 
typical rate of spin-injection into the node is faster than 
the spin-flip relaxation rate, so that the right hand side of 
(|l|) can be set to zero. 

The current through the normal metal node is invari- 
ant with respect to a flip in the magnetization direction 
of ferromagnet Fl or F2: I{9) = I{d + vr). However the 
chemical potential of ferromagnet F2 changes during the 
same process since it is sensitive to the magnitude and 
direction of the spin-accumulation on the normal metal 
node: fi2{d) ^ fJ,2{S + tt). A spin-'resistance' Rs{d) can be 
defined as the ratio between the difference in the chemi- 
cal potential of ferromagnetic F2 when ferromagnet Fl or 
ferromagnet F2 is flipped: Rs = (F2(6'-f vr) - V2{9))/I{9). 
In the coUinear configuration (6 = 0) the spin-resistance 
Rs{0 = 0) is thus the ratio between the difference in the 
chemical potential of ferromagnetic F2 when its magneti- 
zation is parallel (/if) and anti-parallel (/i^^) to the mag- 
netization of ferromagnet Fl and a current (/) passes from 
ferromagnet Fl to the normal metal reservoir N. With the 
aid of the general conductances, valid for arbitrary con- 
tacts, we solve for the non-equilibrium distribution func- 
tion on the normal metal node (|l|) under the condition 
that no particle current enters ferromagnet F2. Using the 
solution for the non-equilibrium distribution function we 
find the current (^) through the system and subsequently 
the non-equilibrium chemical potential of ferromagnet F2. 

Let us first discuss the results in the collinear config- 
uration, = and TT. The spin-resistance can be simply 
expressed as 



Rs{0 = 0) = 



2P1P2 



Gi(l-P2) + G3Ar + G2(l-P|) 



, (47) 



and is independent of the relative mixing conductances rji 
and 772 that are only relevant for the transport proper- 
ties in systems with non-collinear magnetization configu- 
rations. The spin-resistance is proportional to the product 
of the polarizations of the contacts to ferromagnet Fl and 
ferromagnet F2. In order to measure a large effect of the 
spin- accumulation, e.g. a large spin-resistance, highly re- 
sistive contacts should be used. On the other hand, the re- 
sistance has to be small enough so that the transport dwell 
time is shorter than the spin-flip relaxation. The simple re- 
sult (^^ covers a large class of experiments, since we have 
not specified any details about the contacts between the 
reservoirs and the normal metal node. It is noted, though. 



that Eq. (|4^) is only valid for a normal metal node that 
is smaller than the spin-diffusion length and can therefore 
not be applied directly to Johnson's experiment Q. 

We can understand that the present results are quite 
different from those of Ref. |2^ as follows. The general 
formulation in terms of transmission probabilities of Geux 
et al. is exact. However, in order to include the effects of 
spin-relaxation the transmission probabilities were treated 
as pair-wise resistors between the reservoirs. This corre- 
sponds to an equivalent circuit in which resistors con- 
nect the three reservoirs in a "ring" topology. The present 
model, on the other hand, can be described by a "star" 
configuration circuit, in which all resistors point from the 
reservoirs to a single node. The present model is more 
accurate when the contacts dominate the transport prop- 
erties, whereas Geux's model is preferable when the resis- 
tance of the normal metal island is important. Effectively, 
Johnson's thin fllm device appears to be closer to the star 
conflguration. 

Let us now proceed to discuss the results when the 
magnetization directions are non-collinear. The analytical 
expression for the spin-resistance is much simpler when 
the two contacts Fl-N and F2-N are identical, Gi = G2 = 
G, Pi = P2 = P and 771 = 772 = 77. Furthermore we 
disregard the imaginary part of the mixing conductance 
which is very small or zero in the model calculations of 
tunnel, ballistic and diffusive contacts presented in this 
paper as well as in recent first-principle band-structure 
calculations!^. The spin- resistance then has the simple 
form 



Rs 



2(G3jv + 2Gr))P^ cos(6i) 



{GsN + G(l + 7/ - P2))2 _ Q2 cos2(0)(l - 77 - p2)2 

The spin-resistance is an even function of the relative an- 
gle between the magnetization directions 6 and we recover 
the result (^^ when 9 = 0. The spin-resistance vanishes 
when the magnetizations are perpendicular 9 = 7r/2 as 
expected from the symmetry of the systems. The angu- 
lar dependence is approximately proportional to cos(6') 
when the relative mixing conductance is not too large, 
77 « 1. For larger mixing conductances 77 ^ 1 the spin- 
accumulation on the normal metal island is strongly sup- 
pressed in the perpendicular configuration 9 = ^/2 due 
to the large transport rates for spins between the nor- 
mal metal node and ferromagnet F2. Consequently the 
spin-resistance is small and only weakly dependent on the 
relative angle around 9 = tt/2. 

Another novel three-terminal device comprising of three 
ferromagnetic reservoirs (the "spin-flip transistor" ) , which 
utilizes the added functionality provided by non-collinear 
magnetization directions, was introduced in Ref. ||Tq]. 



5 Conclusion 

We developed a mesoscopic circuit theory of spin-transport 
in multi-terminal hybrid ferromagnet - normal metal sys- 
tems starting from microscopic principles. Based on con- 
servation of spin-current on each node the circuit theory 
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is parameterized by the conductances of the contacts, viz. 
two spin-dependent conductances and a (complex) mix- 
ing conductance. The latter is a novel concept relevant 
for transport in systems with non-coUinear magnetization 
configurations. Explicit expressions for the conductances 
for diffusive, ballistic and tunnel contacts have been de- 
rived. The circuit theory leads to simple and quite gen- 
eral results for the conductance of two-terminal and three- 
terminal devices, like Johnson's spin-transistor. For two- 
terminal systems a universal lower limit for the current 
modulation as a function of the relative magnetization 
has been found. 

After completion of this work a complementary ap- 
proach to spin-transport in ferromagnetic-normal metal 
systems starting from the scattering matrices was pre- 
sented in Ref . p4| . Using random matrices to describe the 
scattering within disordered normal metal nodes equiva- 
lent results to our (0) were obtained. 
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A Transformation from a non-collinear to a 
collinear configuration 

We consider the transmission and reflection matrices be- 
tween a normal metal and a ferromagnet. The Schrodinger 
equation is 

H{v),p{v) = £;^(r), 

where ipir) is a two-component spinor. The Hamiltonian 
is 



--^V^i-f y,(r)a, +P;(r) 
2 m 



(49) 



We consider transport between a normal metal and a uni- 
form ferromagnet, so that the magnetization direction m 
is a spatially independent unit vector. In ( |4S| ) T4(r) de- 
notes the spin-dependent potential and K(r) is the scat- 
tering potential of the contact. The direction of the mag- 
netization is represented by the angles 9 and 95 as m = (sin 6 
The hermitian and unitary matrix \J that diagonalizes the 
spin-dependent potential is 



cos(0/2) sin(6'/2)e-*'^/2 
sin(6l/2)e*'^/2 _cos(6l/2) 



The spin-dependent potential vanishes in the normal metal, 
V^s(r) = for X < xi and attains a constant value in the 



ferromagnet T4(r) = Vs for x > Xr- The contact is repre- 
sented by the scattering potential 



v;(r) = 



^sf(r) 
^Jf(r) n(r) 



(50) 



and attains the bulk values within the normal metal and 
the ferromagnet for x < xi and x > a;,., respectively. The 
off-diagonal terms in ( po|) represent the exchange poten- 
tials due to a non-collinear magnetization in the contact, 
spin-orbit interaction or spin-flip scatterers. The Hamilto- 
nian ( ^Of ) can be diagonalized in spin-space by 

^P{r) = U(j){r) . (51) 

The Schrodinger equation for the spinor (/)(r) is 

1 



2m 



V'l + Vsir)a, + Vc{r)-E 



0. 



Let us now consider an incoming wave from the normal 
metal in the transverse mode n and with spin s collinear to 
the magnetization in the ferromagnet. The wave function 
in the normal metal is 



XNip) 



If ik"x I mn 



-ik'^x-\ 



(52) 



where ^| = (1, 0) and — (0, 1) are the spinors, X^{p) is 
the transverse wave function, fc"* is the longitudinal wave 
vector for mode m and is the reflection matrix from 
state ns to state ms'. We would like to transform the result 
for the reflection matrix into a basis with arbitrary spin 
quantization axis. To this end we introduce the incoming 
spinor wave function 



C(r) 



' , /l.m 



(53) 



Using the transformation ( pi] ) we can also write the wave 
function spinor in terms oiMthe basis states 0(r) as 



C(r) = f7VC(r)a.., 



(54) 



where a^s are expansion coefficients to be determined by 
equating ([5^) and (Q). We thus find that 



cos 95, sm a sm 1^9, cos ( 
and 



7 , ^s'a'l Ca'a 



(55) 



Disregarding spin-flip processes in the contact the trans- 
formation of the reflection matrix can be written as 
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where the spin-projection matrices ui and U| are defined 
m (|) and (|). Spin-flip processes can also be included by 
using the general transformation (|55|), but the reflection 
matrix f"™ can then not be expressed in terms of the 
spin-projection matrices only. 

We can perform a similar calculation in order to find 
the transformation of the transmission coefficients from 
the ferromagnet into the normal metal. In the basis where 
the spin-quantization axis is coUinear with the magneti- 
zation the incoming wave from the ferromagnet is 



X'Psip) 



(56) 



where x'rnip) is the (spin-dependent) transverse wave func- 
tion and fc™ is the spin-dependent Fermi wave- vector. The 
outgoing wave into the normal metal is 



By transforming the outgoing wave into an arbitrary mag- 
netization direction according to (|5l|), we see that the 
transmission coefficient from a state with spin s collinear 
to the magnetization direction in the ferromagnet to a 
state with spin s' collinear to the spin-quantization axis 
along the z-direction is 



TTj.mn 



In the absence of spin-flip scattering in the contact, we 
have 

The current in the normal metal is (for a given energy 
shell) (note that we associate the first index with 'f' and 
the second index with !Z't) 



-/ = Mf 

e 



N 



The contribution from the transmission probability to the 
spin-current is therefore 



tF ^ ,lmn rF 



fnm \ t 
f3S ) 



nm/3 



nmf3 



fmn I 
I 



Whereas the contribution from the transmission probabil- 
ity becomes 



where the spin-dependent conductance is 

4./mn 1 2 



nm 



Similarly, the contribution from the normal metal is 

ss' 



where we have used the unitarity of the scattering matrix 

so that M - E„™ = Enm 1^'"™!' 

conductance is introduced as 



and the mixing 



G 



n 



e 



mn \ nm 
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Fig. 1. A many-terminal circuit consisting of normal metals 
(N) and ferromagnets (F) with arbitrary magnetization direc- 
tions. The normal metals (ferromagnets) are coupled to adja- 
cent nodes by contacts which determine the resistance of the 
system. 

Fig. 2. A contact between a ferromagnetic node and a normal 
metal node. The current is evaluated at the normal metal side 
(dotted line). The transmission coefficient from the ferromag- 
net to the normal metal is t' and the reflection matrix from 
the normal metal to the normal metal is r. 

Fig. 3. A diffusive contact between a ferromagnetic node and 
a normal metal node. The length of the ferromagnetic (normal 
metal) part of the contact is {L'^), and the conductivity of 
the ferromagnetic (normal metal) is (o"'^)- The cross-section 
of the contact is A. 

Fig. 4. The two-terminal device comprising a normal metal 
node attached to two ferromagnetic reservoirs (Fl and F2) 
with arbitrary relative magnetization direction 9. A source- 
drain bias V is applied between the ferromagnetic reservoirs 
and a current I flows between the two reservoirs. The contact 
between ferromagnetic Fl (F2) and the normal metal node is 
characterized by the conductances G}, G{, and Gl^ {Gl, Gj, 
and Gl^). 



Fig. 6. The three terminal Johnson spin-transistor. A bias 
voltage is applied between a ferromagnet and a normal metal 
and a current flows between the same reservoirs through a 
normal metal node. The potential on another ferromagnet is 
measured when its magnetizations is parallel or anti-parallel to 
the first ferromagnet. The contact between the first (second) 
ferromagnet and the normal metal node is characterized by the 
total conductance Gi (G2), and the polarization Pi (P2). The 
conductance between the normal metal node and the normal 
metal reservoir is Gsat. 
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Fig. 5. The current through the two-terminal device as a func- 
tion of the relative angle 9 between the magnetization direc- 
tions in Fl and F2. The current is normalized by the current 

p 

in the parallel configuration I — GV/2 and reaches the min- 
imum J = (1 - PV(1 + 9sf))GV/2 = (1 - P^g)GV/2 when the 
magnetizations are anti-parallel {9 — tt). The lowest solid line 
is for the minimum relative mixing conductance \ri = 1|. This 
line forms the lowest possible universal curve for the conduc- 
tance vs. relative angle 9, since all other values of the mixing 
conductance lies above this line. 
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